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X-RAY  EMISSION  FROM  CHj  LASER-PRODUCED 
PLASMAS  AT  101*  W/cm2 


INTRODUCTION 

Experimental  X-ray  spectra  have  been  used  to  demonstrate  the  existence  of  suprathermal 
electrons1,  ion  acoustic  turbulence2  and  self-generated  magnetic  fields3.  The  first  two 
phenomena  imply  nonthermal  electron  distribution  functions  whereas  the  third  one  results  in 
the  production  of  energetic  thermal  electrons.  Davis  and  Rogerson4  have  shown  that  it  is 
quite  difficult  to  distinguish  between  these  different  types  of  electron  distribution  functions 
from  Bremsstrahlung  measurements,  the  main  difference  being  at  high  photon  energies  where 
statistical  measurement  errors  are  usually  large. 

Since  the  experimental  data  does  not  exclude  the  magnetic  field  interpretation,  we  look 
further'  into  it.  We  present  more  details  of  the  physical  model  and  show  where  the  different 
components  of  the  X-ray  spectra  originate.  We  illustrate  the  dependence  of  the  computed 
spectra  on  different  physical  parameters.  We  present  the  dynamics  a two-dimensional  numeri- 
cal model  predicts  for  a laser-slab  target  interaction.  We  finally  include  tests  of  X-ray  spectra 
sensitivity  to  parameters  such  as  grid  size,  flux-limiting  factors  and  different  laser  absorption 
mechanisms. 
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MODEL 

The  basic  model  has  been  described  previously  3.  Some  of  its  unique  features  are  that  it 
uses  an  equation  of  state 5 rather  than  an  ideal  gas  law  to  describe  the  plasma  produced,  it  pro- 
vides a complete  treatment  of  magnetic  fleld  effects,  including  the  tensor  thermal  conductivity. 
Also,  the  2-dimensional  code  does  not  contain  any  adjustable  parameters. 

In  addition  to  the  features  mentioned  in  Ref.  13],  a simple  form  for  radiation  pressure6 
has  been  included  in  the  momentum  equation  and  a more  complete  thermal  conductivity  ten- 
sor has  been  included.  In  particular,  the  Peltier  heat  conduction  term  q — /3  • u has  been  ad- 
ded where  u is  the  electron  drift  velocity  and  /3  is  given  by 


where  the  /3*s  and  x — u>t  are  defined  in  ref.  7.  The  cross-conductivity  in  the  Fourier  heat 


conduction  term 


has  also  been  included.  We  wish  to  emphasize  the  fact  that  we  use  the  classical  conductivity 


and  do  not  flux-limit  our  transport.  As  we  show  below,  results  with  and  without  flux-limiting 


are  nearly  identical.  This  may  be  interpreted  as  a demonstration  that  our  fluid  model  is  self- 


consistent.  Should  the  flux-limited  and  non-flux-limited  results  be  significantly  different,  then 


the  hypotheses  on  which  the  fluid  model  was  derived  would  be  violated. 


The  equation  for  the  magnetic  fleld  as  reported  in  this  work  does  not  contain  the  thermal 


force  whose  effects  are  reported  in  detail  in  ref.  8.  The  spectral  results  reported  here  are  ob 


tained  with  free-free  X-ray  diagnostics  which  have  been  added  to  the  code  by  breaking  the 


NRL  MEMORANDUM  REPORT  363S 


bremsstrahlung  radiation  rate  into  17  frequency  intervals: 


Pbrtm  (£\  • E2  ^ 


p -< i 

rbl  e 


Here  Pb,  and  T,  are  the  local  total  bremsstrahlung  radiation  rate  and  temperature  obtained 
from  the  local  ion  and  energy  densities  through  the  equation  of  state.  The  frequency  intervals 
were  chosen  to  correspond  to  the  range  of  experimental  interest  and  the  radiation  output  in 
each  frequency  range  was  integrated  over  space  and  time  to  obtain  the  total  output.  For  the 
X-ray  energies  considered,  it  was  assumed  that  half  the  radiation  escaped  from  the  plasma  into 
2n  steradians,  and  the  remaining  half  was  deposited  in  the  solid  region. 


STANDARD  CASE 


All  of  the  results  which  are  presented  in  this  section  correspond  to  the  following  input 
parameters:  1.2  J of  1.06  n light  deposited  in  a 12  m radial  focal  spot  in  a pulse  with  a FWHM 
of  20  psec.  The  pulse  is  Gaussian  in  space  and  time  and  has  a peak  flux  of  1.2  10 16  W/cm2. 
No  prepulse  precedes  the  main  pulse.  The  target  is  a plane  polyethylene  slab  • and  a steep  gra- 
dient density  has  been  set  up  in  front  of  it  in  order  to  start  the  calculations.  Figure  1 depicts 
the  initial  laser  and  target  conditions.  The  laser  parameters  are  typical  of  the  experiments  per- 


formed with  the  NRL  laser  facility3. 

In  the  following,  we  present  in  detail  the  numerical  results  obtained  with  the  model 
described  in  the  previous  section  on  a 12  m by  12  m grid-  Then,  we  investigate  different  physi- 
cal effects  such  as  variation  of  the  initial  density  gradient,  the  ion  atomic  mass,  the  absorption 
coefficient  and  of  the  inclusion  of  a flux-limiter.  Finally  we  present  partial  results  obtained  for 
the  same  physical  case  with  a by  V grid. 

Figure  2 thematically  pictures  the  target  density  and  temperature  profiles  and  illustrates 
some  features  of  these  profiles  which  will  be  used  to  explain  the  magnetic  field  dynamics 
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below.  Figures  3 through  S are  profiles  of  fluid  variables  along  the  axis  for  selected  times  near 
peak  intensity  and  later  in  the  pulee  (from  2 to  72  peec  after  the  peak  of  the  laser  pulse).  In 
these  figures,  the  maximum  in  temperature  occurs  near  the  critical  surface  (where  the  absorp- 
tion is  maximum)  and  it  follows  rather  closely  the  laser  power  variation  with  time.  It  reaches 
M keV  2 peec  after  the  peak  of  the  laser  pulse  and  drops  to  about  1 keV  72  psec  after  that 
peak.  Note  that  at  i — 21  psec  after  the  pulse  maximum,  a series  of  maxima  appears  in  the 
temperature  profile.  These  maxima  might  be  an  indication  of  an  acoustic  instability  driven  by 
the  magnetic  source  terms. 

The  density  profile  has  been  reshaped  by  the  laser  pulse  by  the  time  of  the  peak  of  the 
laser  pulse  and  shows  a very  large  gradient  in  the  vicinity  of  the  critical  surface.  The  density 
jumps  by  a factor  of  50  to  100  over  two  cells  or  24p.  This  sharp  gradient  is  moving  inwards 
into  the  target  whereas  a rarefaction  shock  propagates  outwards,  towards  the  laser  beam.  This 
feature  might  suggest  that  a steady-state  density  gradient  forms  rapidly  near  the  critical  surface 
and  that  the  experimental  results  are  thus  very  weakly  dependent  on  the  presence  of  a 
prepulse.  This  feature  also  explains  the  weak  dependence  of  the  numerical  results  with  the 
values  of  the  intitial  density  gradient.  We  remark  here  that  special  numerical  techniques9  are 
used  here  which  accurately  preserve  the  physics  of  such  shock  processes 10. 

The  self-generated  B-field  is  also  shown  in  those  figures  (off-axis  at  r - 12  m)  and  we  see 
that  its  maximum  value  (when  the  thermal  force  is  not  included)  remains  almost  constant  at 
over  2 MG  for  the  first  18  psec.  This  magnetic  field  is  a dominant  feature  of  the  interaction 
region  and  once  generated,  difftises  both  into  the  target  and  towards  the  laser.  The  maximum 
B-field  occurs  near  the  critical  surface  and  field  reversals  increase  with  time.  We  recall  that 
these  field  reversals  diminish  with  the  addition  of  the  thermal  force  term  *. 


The  complicated  B-fleld  profile  results  from  the  structure  of  the  density  and  temperature 
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gradients.  Note  however  that  the  dominant  term  is  dn/dz  8T/dr.  The  positive  radial  density 
gradient  term  due  to  the  presence  of  the  initial  density  gradient  gives  a much  smaller  contribu- 
tion to  the  field  source  term  (dn/dr  8T  lb z)  and  thus  shows  that  the  initial  density  gradients 
only  perturbs  slightly  the  solution.  The  maximum  value  of  the  B-fteld  remains  at  a high  level 
during  the  laser  pulse,  because  as  we  have  seen,  dn/dz  is  approximately  constant  during  the 
laser  interaction  as  well  as  dT/dr  which  remains  also  constant  due  to  the  reduced  conductivity. 
It  is  also  apparent  from  Fig.  5 that  the  B-field  is  carried  off  with  the  plasma  as  the  gradients  in 
the  expanding  shock  move  away. 

The  parameter  <u„  rei  is  a measure  of  the  reduction  in  the  perpendicular  thermal  conduc- 
tion. Its  peak  follows  the  maximum  in  magnetic  field  and  reaches  a value  greater  than  10 J. 
The  consequent  reduction  in  heat  conduction  at  that  time  varies  from  a factor  of  10 J for  the 
cross-conductivity  to  10 6 for  the  perpendicular  one.  It  is  this  reduction  in  heat  conduction 
which  acounts  for  the  high  temperatures  observed  numerically  and  for  the  hard  X-rays  which 
are  also  observed  experimentally. 

The  total  radiation  power  output  peaks  just  behind  the  critical  surface  where  the  density 
is  large  and  the  temperature  still  above  a few  hundreds  of  eV.  At  early  times,  a strong  emit- 
ting region  is  also  seen  associated  with  the  outward  moving  shock,  but  as  the  shock  wave 
moves  to  lower  densities  and  cools  by  expansion,  this  strongly  emitting  region  disappears.  The 
maximum  radiated  power  region  is  shown  as  a function  of  time  in  Fig.  6.  At  first  it  occurs  in- 
side the  spike  the  laser  beam  forms  inside  the  target.  Then  it  spreads  outwards  to  reappear  on 
the  axis.  This  pattern  can  be  explained  in  terms  of  diffusion  and  fluid  motion.  The  heat  is  car- 
ried radially  by  convection  and  conduction,  taking  the  emitting  region  with  it  but  when  the 
density  drops  too  low,  the  maximum  emitting  region  reappears  on  the  axis. 

Fig.  7 to  12  show  temperature  and  velocity  contours  for  various  times.  Velocities  are 
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directed  at  first  toward  the  laser  beam  and  then  spread  out  radially.  Also  the  front  of  the  ex- 
pansion accelerates  from  4 x 10 7 cm/sec  at  / - 18  pscc  to  10 8 cm/sec  at  t - 88.5  psec  which 
correspond  to  average  ion  energies  going  from  7.7  keV  to  48.4  keV. 

Before  interpreting  the  spatially  time-integrated  X-ray  spectra,  we  have  plotted  the  X-ray 
power  emitted  between  1 and  5 keV  and  above  5 keV  in  Fig.  13.  We  see  that  the  maximum  of 
the  harder  X-rays  comes  earlier,  around  30  psec  after  the  beginning  of  the  interaction  when 
the  laser  pulse  is  still  on.  This  indicates  clearly  that  the  hard  X-rays  are  due  to  the  high  tem- 
peratures generated  in  the  plasma  near  the  peak  of  the  laser  pulse.  Lower  energy  X-rays  come 
later  when  the  temperature  drops  due  to  the  decrease  in  laser  power  and  to  heat  diffusion.  Re- 
cently, these  time-resolved  X ray  spectra  have  just  been  obtained  experimentally 1 1 and  should 
provide  the  basis  for  a more  detailed  comparison  than  that  one  provided  below. 

The  spatially  and  time  integrated  X-ray  spectrum  is  shown  in  Fig.  14  along  with  the  ex- 
perimental spectrum  obtained  at  NRL^.  By  comparison,  a spectrum  obtained  with  the  magnet- 
ic field  artificially  set  to  0 is  also  shown.  The  hard  X-rays  have  been  already  explained  in 
terms  of  the  high  temperatures  (tens  of  keV)  generated  during  the  laser  interaction.  However, 
these  calculated  intensities  are  larger  than  the  measured  ones.  Although  the  inclusion  of  the 
thermal  force  terms  reduce  the  X-ray  intensities,  other  reasons  can  be  given  to  explain  the 
larger  X-ray  intensities  found  numerically.  We  shall  mention  two  reasons  here.  (1)  Electrons  of 
100  keV  for  example  have  a Larmor  radius  of  8 m in  a 1 MG  B-field,  which  is  comparable  with 
the  size  of  the  high  B-field  region.  Thus,  the  electrons  would  see  an  inhomogeneous  magnetic 
field  and  drift  out  of  the  hot  region.  (2)  Another  fact  is  that  for  T,  -40  keV  and  ne  - 10 21 
cm  ~3,  the  electron  self-collision  time  is  400  psec  which  is  much  larger  than  the  characteristic 
times  in  the  experiment.  So,  although  electrons  might  be  produced  at  40  keV,  they  do  not 
maxwellize  on  the  time  scales  of  the  experiment  and  as  a result,  they  contribute  to  an  overes- 
timated hard  X-ray  contribution.  This  very  hard  part  of  the  spectrum  is  also  that  one  where 


NRL  MEMORANDUM  REPORT  3638 


the  experimental  unceriainty  is  the  largest  since  the  number  of  photons  in  this  range  is  orders 
of  magnitude  below  than  in  lower  energy  ranges. 

» i . 

The  soft  component  of  our  calculated  X-ray  spectrum  is  much  more  difficult  to  compare 

I I , . 

accurately  with  experiment  because  sharp  gradients  in  both  densities  and  temperature  occur 
over  1 or  2 cells.  The  resolution  is  poor  for  that  region  of  the  laser  target  interaction  and  it  will 
be  seen  later  how  this  part  of  the  spectrum  varies  with  the  grid  size.  In  any  case,  the  agree- 
ment is  reasonably  good  for  energies  above  1 keV  which  correspond  to  temperatures  just  on 
the  border  of  that  sharp  gradient  region.  One  feature  of  the  experimental  X-ray  spectrum  does 
not  show  in  the  numerical  curve;  it  is  the  dip  in  the  spectrum  between  2 and  10  keV.  This  dip 
seems  to  indicate  a decoupling  mechanism  between  low  energy  electrons  and  higher  energy 
electrons.  In  particular,  the  higher  energy  electrons  might  not  be  suited  for  a temperature 
representation  (since  they  do  not  have  time  to  maxwellize)  and  they  would  not  obey  either  a 
diffusion  equation.  These  remarks  do  not  reflect  on  the  origin  of  suprathermal  electrons  but 
suggest  that  once  suprathermal  have  been  identified  in  the  physical  interaction,  their  transport 
might  have  to  be  done  differently. 

When  the  magnetic  fields  are  artificially  set  to  0,  the  hard  X-ray  component  of  the  spec- 

< f • '•*  • - * 

trum  disappears  almost  completely  since  for  h v — 40  keV,  the  X-ray  intensity  is  4 orders  of 
magnitude  smaller.  This  constitutes  strong  evidence  for  the  role  played  by  large  B-fields  in 
producing  the  high  temperatures  needed  to  explain  the  hard  X-ray  experimental  data. 

RESULTS 

• • , ‘ : . A ’ *-  , 

We  now  discuss  specific  conclusions  obtained  from  comparisons  of  calculations  with  this 
model,  and  their  dependences  on  physically  interesting  parameters.  We  consider  the  influence 
of  the  initial  density  gradient,  the  ion  mass,  the  laser  absorption  mechanism,  the  importance  of 

fine  structure,  and  the  presence  of  flux-limiting. 
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As  mentioned  previously,  the  initial  density  gradient  is  not  playing  a significant  role  in 
the  value  of  the  integrated  X-ray  spectrum.  Results  shown  in  Fig.  IS  with  density  gradient 
length  varying  from  20/i  to  41m  illustrates  this  result  since  the  variation  in  X-ray  intensity  is  a 
factor  of  2.S  at  most.  The  above  variations  in  density  gradient  lengths  correspond  to  a dip  in 
densities  varying  by  4 orders  of  magnitude  over  the  distance  comprised  between  the  boundary 
and  the  critical  surface.  This  result  reinforces  the  conclusion  that  the  laser  light  produces  its 
own  density  profile  in  the  interaction  region  and  that  the  X-ray  spectrum  is  insensitive  to  the 
presence  of  a prepulse.  The  density  structure  near  the  critical  surface  might  also  explain  some 
features  of  backscattered  light. 

The  model  which  we  have  used  to  describe  the  laser  target  interaction  is  a one  fluid 
model  whereas  two  ion  species  make  up  the  target  material:  carbon  and  hydrogen.  Leaving 
aside  the  question  of  ion  separation  and  acceleration,  computations  were  made  with  carbon  or 
polyethylene  mass  for  the  ions.  The  resulting  spatially-and  time-integrated  X-ray  spectra  are 
shown  in  Fig.  16  These  two  curves  are  very  similar  since  a factor  of  at  most  3 separates  the 
X-ray  intensities  over  the  entire  X-ray  energy  range.  In  view  of  the  similarity  of  these  results 
and  because  the  polyethylene  case  was  much  longer  to  run  than  the  carbon  case  (since  the 
velocities  are  larger  for  lighter  ions  and  as  a consequence,  the  time  step  is  much  smaller),  all 
other  cases  shown  were  run  with  the  carbon  ion  mass. 

Absorption  by  an  evanescent  wave  in  the  overdense  region  was  also  considered  as  an  al- 
ternate energy  deposition  mechanism.  An  arbitrary  fraction  (50%)  of  the  energy  reaching  the 
critical  surface  was  assumed  to  go  into  the  evanescent  wave.  Since  more  energy  gets  behind 
the  critical  surface  and  less  energy  heats  the  critical  surface,  the  maximum  temperature  is  ex- 
pected to  be  lower  than  for  the  standard  case  but  higher  in  the  overdense  region.  This  should 
lead  to  a loss  of  intensity  for  hard  X-rays  and  increase  in  intensity  for  soft  X-rays.  Figure  17 

shows  these  trends  where  variations  are  at  most  a factor  of  4. 
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The  most  severe  limitation  on  the  results  presented  here  is  undoubtedly  the  spatial  reso- 
lution in  the  vicinity  of  the  critical  surface  where  large  gradients  occur.  A limited  attempt  was 
made  at  assessing  the  influence  of  the  grid  size  on  the  solution.  A run  was  performed  up  to 
the  time  of  the  peak  of  the  laser  pulse  with  a 4/*x  4 n grid  instead  of  the  12  n x 12  n grid 
used  in  the  standard  case.  Profiles  for  density,  temperature,  B-field  show  more  structure  at  the 
same  time  and  the  integrated  X-ray  spectra  is  shown  in  Fig.  18.  This  spectrum  shows  a 
difference  of  a factor  of  2 in  the  hard  X-ray  region  and  an  increase  by  a factor  of  2 in  the  soft 
X-ray  range.  This  increase  confirms  the  need  for  a more  detailed  solution  in  regions  of  large 
steep  gradients. 

In  conclusion,  although  the  details  of  the  interaction  might  change  slightly  due  to  the 
spatial  resolution,  the  integrated  X-ray  spectra  which  reflect  gross  features  like  maximum  tem- 
perature obtained  and  extent  of  heated  regions  agree  reasonably  well  and  do  not  raise  questions 
about  different  phenomena  taking  place  as  a result  of  the  interaction. 

Because  of  large  gradients  in  temperature,  the  thermal  flux  may  become  larger  than  1/2 
nt  me  v,3*.  Runs  with  and  without  flux  limiter  were  performed  and  results  appear  in  Figs.  19 
and  20.  In  Fig.  19,  electron  temperature  profile  on  the  axis  at  the  peak  of  the  laser  pulse  is 
shown  for  the  cases  B — 0 and  B ^ 0,  with  and  without  flux  limiter.  It  is  seen  that  the 
difference  between  the  cases  with  and  without  flux-limiter  is  relatively  small.  The  correspond- 
ing X-ray  spectra  are  shown  in  the  following  figure.  The  higher  intensities  for  the  hard  X-ray 
follow  from  the  higher  temperatures  observed  in  the  previous  figure.  The  effect  of  the  flux- 
limiter  is  very  small  on  our  results  and  indicates  that  the  calculations  are  performed  well  within 
the  limits  of  fluid  theory. 
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CONCLUSION 

We  have  found  good  agreement  between  numerical  and  experimental  X-ray  spectra  for 
experimentally-measured  plasma  parameters,  during  the  interaction  when  a CH2  target  is  irra- 
diated by  1016W/cm2  of  1.06  light.  We  have  shown  that  this  agreement  is  preserved  when 
initial  conditions,  ion  mass,  laser  energy  deposition  and  grid  size  are  changed  but  not  when  the 
self-generated  B-fields  inside  the  laser-produced  plasma  was  artificially  set  to  0.  These  results, 
supplemented  by  results  in  ref.  8,  indicate  that  B-fields  which  have  been  measured  directly  and 
indirectly  represent  one  of  the  most  probable  explanations  for  the  observation  of  the  hard  X- 
ray  spectrum. 

We  have  shown  also  that  the  hard  X-ray  pulse  follows  very  closely  the  laser  pulse  in 
time,  a fact  which  could  be  checked  eventually  in  an  experiment.  We  have  also  shown  that  a 
more  accurate  prediction  of  soft  X-rays  requires  improved  resolution  in  the  region  of  steep  gra- 
dients. We  draw  attention  to  the  small  influence  played  by  a flux-limiter  in  our  calculations. 
Finally,  we  observe  that  oscillations  in  the  temperature  profiles  could  come  from  an  acoustic 
instability  driven  by  the  magnetic  source  terms. 
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Schematic  of  the  radiation  contour  intensities  at  various  limes  during  the  interaction. 
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RESULT  WITHOUT  FLUX-UMTER 
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Numerical,  experimental  and  magnetic  field-less  X-ray 
spectra  as  a function  of  photon  energy. 
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Fig.  I S — Variition  or  X-ray  spectra  at  t — 300  psec  as  a function  of  the  ratio 
%/Ami  Fit-  I)  — 'Standard'  case  corresponds  to  K - 10~V 


